Linear-scaling first-principles molecular dynamics of complex biological systems with the CONQUEST code

The recent progress of linear-scaling or $\mathcal{O}(N)$ methods in density functional theory (DFT) is remarkable. In this paper, we show that all-atom molecular dynamics simulations of complex biological systems based on DFT are now possible using our linear-scaling DFT code Conquest. We first overview the calculation methods used in Conquest and explain the method introduced recently to realise efficient and robust first-principles molecular dynamics (FPMD) with $\mathcal{O}(N)$ DFT. Then, we show that we can perform reliable all-atom FPMD simulations of a hydrated DNA model containing about 3400 atoms. We also report that the velocity scaling method is both reliable and useful for controlling the temperature of the FPMD simulation of this system. From these results, we conclude that reliable FPMD simulations of complex biological systems are now possible with Conquest

The PCI Interface for GRAPE Systems: PCI-HIB

We developed a PCI interface for GRAPE systems. GRAPE(GRAvity piPE) is a special-purpose computer for gravitational N-body simulations. A GRAPE system consists of GRAPE processor boards and a host computer. GRAPE processors perform the calculation of gravitational forces between particles. The host computer performs the rest of calculations. The newest of GRAPE machines, the GRAPE-4, achieved the peak performance of 1.08 Tflops. The GRAPE-4 system uses TURBOChannel for the interface to the host, which limits the selection of the host computer. The TURBOChannel bus is not supported by any of recent workstations. We developed a new host interface board which adopts the PCI bus instead of the TURBOChannel. PCI is an I/O bus standard developed by Intel. It has fairly high peak transfer speed, and is available on wide range of computers, from PCs to supercomputers. Thus, the new interface allows us to connect GRAPE-4 to a wide variety of host computers. In test runs with a Barnes-Hut treecode, we found that the performance of new system with PCI interface is 40% better than that of the original system.Comment: 15 pages, 10 Postscript figures, 3 tables, Latex, submitted to Publications of the Astronomical Society of Japan. corrected figure 2 which contained non standard font

Interfacial Layer Breaker: A Violation of Stokes' Law in High-Speed Atomic Force Microscope Flows.

Structured water near surfaces is important in nonclassical crystallization, biomineralization, and restructuring of cellular membranes. In addition to equilibrium structures, studied by atomic force microscopy (AFM), high-speed AFM (H-S AFM) can now detect piconewton forces in microseconds. With increasing speeds and decreasing tip diameters, there is a danger that continuum water models will not hold, and molecular dynamic (MD) simulations would be needed for accurate predictions. MD simulations, however, can only evolve over tens of nanoseconds due to memory and computational efficiency/speed limitations, so new methods are needed to bridge the gap. Here, we report a hybrid, multiscale simulation method, which can bridge the size and time scale gaps to existing experiments. Structured water is studied between a moving silica AFM colloidal tip and a cleaved mica surface. The computational domain includes 1,472,766 atoms. To mimic the effect of long-range hydrodynamic forces occurring in water, when moving the AFM tip at speeds from 5 × 10-7 to 30 m/s, a hybrid multiscale method with local atomistic resolution is used, which serves as an effective open-domain boundary condition. The multiscale simulation is thus equivalent to using a macroscopically large computational domain with equilibrium boundary conditions. Quantification of the drag force shows the breaking of continuum behavior. Nonmonotonic dependence on both the tip speed and distance from the surface implies breaking of the hydration layer around the moving tip at time scales smaller than water cluster formation and strong water compressibility effects at the highest speeds

Coupled Replicator Equations for the Dynamics of Learning in Multiagent Systems

Starting with a group of reinforcement-learning agents we derive coupled replicator equations that describe the dynamics of collective learning in multiagent systems. We show that, although agents model their environment in a self-interested way without sharing knowledge, a game dynamics emerges naturally through environment-mediated interactions. An application to rock-scissors-paper game interactions shows that the collective learning dynamics exhibits a diversity of competitive and cooperative behaviors. These include quasiperiodicity, stable limit cycles, intermittency, and deterministic chaos--behaviors that should be expected in heterogeneous multiagent systems described by the general replicator equations we derive.Comment: 4 pages, 3 figures, http://www.santafe.edu/projects/CompMech/papers/credlmas.html; updated references, corrected typos, changed conten

A rare case of concomitant huge exophytic gastrointestinal stromal tumor of the stomach and Kasabach-Merritt phenomenon

<p>Abstract</p> <p>Background</p> <p>We report an extremely rare case of concomitant huge exophytic GIST of the stomach and Kasabach-Merritt phenomenon (KMP).</p> <p>Case presentation</p> <p>The patient was a 67-year-old man experiencing abdominal distension since September 2006. A physical examination revealed a 25 × 30 cm hard mass that was palpable in the middle and lower left abdomen minimal intrinsic mobility and massive ascites. Since the admitted patient was diagnosed with DIC, surgery could not be performed. The patient received a platelet transfusion and the DIC was treated. Due to this treatment, the platelet count recovered to 7.0 × 10⁴; tumor resection was performed at 16 days after admission. Laparotomy revealed a huge extraluminal tumor arising from the greater curvature of the stomach that measured 25 × 30 cm and had not ruptured into the peritoneal cavity or infiltrated other organs. Partial gastric resection was performed. The resected mass measured 25 × 25 × 20 cm. In cross section, the tumor appeared hard and homogenous with a small polycystic area. Histopathology of the resected specimen showed large spindle cell GIST with >5/50 HPF (high-power field) mitotic activity. The postoperative course was uneventful, and the coagulopathy improved rapidly.</p> <p>Conclusion</p> <p>Since the characteristic of tumor in this case was hypervascularity with bleeding and necrotic lesions, coagulopathy was thought to be caused by the trapping of platelets within a large vasculized tumor mass.</p

Trends in life science grid: from computing grid to knowledge grid

BACKGROUND: Grid computing has great potential to become a standard cyberinfrastructure for life sciences which often require high-performance computing and large data handling which exceeds the computing capacity of a single institution. RESULTS: This survey reviews the latest grid technologies from the viewpoints of computing grid, data grid and knowledge grid. Computing grid technologies have been matured enough to solve high-throughput real-world life scientific problems. Data grid technologies are strong candidates for realizing "resourceome" for bioinformatics. Knowledge grids should be designed not only from sharing explicit knowledge on computers but also from community formulation for sharing tacit knowledge among a community. CONCLUSION: Extending the concept of grid from computing grid to knowledge grid, it is possible to make use of a grid as not only sharable computing resources, but also as time and place in which people work together, create knowledge, and share knowledge and experiences in a community

Systematic analysis of mitochondrial genes associated with hearing loss in the Japanese population: dHPLC reveals a new candidate mutation

<p>Abstract</p> <p>Background</p> <p>Variants of mitochondrial DNA (mtDNA) have been evaluated for their association with hearing loss. Although ethnic background affects the spectrum of mtDNA variants, systematic mutational analysis of mtDNA in Japanese patients with hearing loss has not been reported.</p> <p>Methods</p> <p>Using denaturing high-performance liquid chromatography combined with direct sequencing and cloning-sequencing, Japanese patients with prelingual (N = 54) or postlingual (N = 80) sensorineural hearing loss not having pathogenic mutations of m.1555A > G and m.3243A > G nor <it>GJB2 </it>were subjected to mutational analysis of mtDNA genes (<it>12S rRNA</it>, <it>tRNA</it>^{<it>Leu(UUR)</it>}, <it>tRNA</it>^{<it>Ser(UCN)</it>}, <it>tRNA</it>^<it>Lys</it>, <it>tRNA</it>^<it>His</it>, <it>tRNA</it>^{<it>Ser(AGY)</it>}, and <it>tRNA</it>^<it>Glu</it>).</p> <p>Results</p> <p>We discovered 15 variants in <it>12S rRNA </it>and one homoplasmic m.7501A > G variant in <it>tRNA</it>^{<it>Ser(UCN)</it>}; no variants were detected in the other genes. Two criteria, namely the low frequency in the controls and the high conservation among animals, selected the m.904C > T and the m.1105T > C variants in <it>12S rRNA </it>as candidate pathogenic mutations. Alterations in the secondary structures of the two variant transcripts as well as that of m.7501A > G in <it>tRNA</it>^{<it>Ser(UCN) </it>}were predicted.</p> <p>Conclusions</p> <p>The m.904C > T variant was found to be a new candidate mutation associated with hearing loss. The m.1105T > C variant is unlikely to be pathogenic. The pathogenicity of the homoplasmic m.7501T > A variant awaits further study.</p